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ABSTRACT 


Dynamic  instability  of  planing  craft  on  calm  water,   porpoising, 
is  a  phenomenon  which  has  not  been  properly  understood.   Empirical 
relations  are  available  for  predicting  the  regime  of  stability.  The  re- 
lations, when  compared,   lead  to  conflicting  design  requirements  to  in- 
crease stability. 

It  is  therefore  desirable  to  develope  a  theoretical  approach  to 
the  problem  so  that  the  effects  of  beam,   deadrise  angle,   etc.  on 
stability  can  be  studied. 

The  results  of  the  investigation  imply  that  a  decrease  in  dead- 
rise  angle,  a  decrease  in  beam  and  an  increase  in  distance  from  LCG 
to  transom  result  In  an  increase  in  stability.  Changes  in  shaft  angle 
and  vertical  height  of  the  center  of  gravity  and  moment  of  inertia  have 
very  little  effect  on  the  stability  of  a  boat  while  it  is  planing.  However 
further  investigation  is  required  to  verify  these  results. 

In  conjunction  with  this  paper,   a  computer  program  was  written 
which  can  be  used  in  the  design  of  planing  craft  to  predict  boat  attitude, 
wetted  surface  area,   drag  and  effective  hor3e  power.  This  program 
will  be  available  for  use  in  the  XIII  Department  library. 
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NOMENCLATURE 

Standard 
Symbol 


b 
C 
C 
D, 


lb 
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KG 


cp 


CG 

lrn 

L 
c 

S 

mr8 
mX 
N 
T 


G 


P 
A 

€ 
6 
X 


Definition 


Beam 

Lift  coefficient  for  prismatic  surface 

Lift  coefficient  for  zero  deadrise  surface 

Drag  force  (lbs) 

Froude  number 


Uo//     r 


Height  of  center  of  gravity  above  base 
line  (ft) 

Non-dimensionalizing  length  (ft) 

Location  of  center  of  pressure  forward 
of  transom  (ft) 

Location  of  center  of  gravity  forward  of 
transom  (ft) 

Mean  wetted  length  (ft) 

Length  of  wetted  chine  (ft) 

Length  of  wetted  keel  (ft) 
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Normal  force  (lbs) 

Thrust 

Mean  velocity  of  flow  past  bottom 

Angle  of  keel  above  horizontal  (deg) 

Non-dimensional  velocity  (fps) 

Wetted  surface  (fta) 

Average  deadrise  angle  (degrees) 

Displacement  (lbs) 

Shaft  angle  (degrees) 

Pitch  angle  (degrees) 

Ratio  of  mean  wetted  length  to  beam 
Note:  Reference  (10)  uses  this  definition 
while  reference  (11)  uses  its  reciprocal. 
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Standard 
Symbol 

P 

T 

(m  -  Z    ) 
w 

Z 
w 


(Z.  +  mXJ 
q  u 

(Z    +  U0ZO 
q  w 

(2B+UoZw) 
(I     -MO 

y      q 

(Mq  ♦  u0  M^) 
<MQ  +  U0  Mw) 

(Mw  +  mXG) 

M 
w 

Mz 


Definition 


Program 
Symbol 


Mass  density  RHO 

Trim  angle  TAU 

Vertical  force  per  unit  vertical  acceleration  Al 
Vertical  force  per  unit  vertical  velocity  Bl 

Vertical  force  per  unit  vertical  displacement  CI 
Vertical  force  per  unit  angular  acceleration  Dl 
Vertical  force  per  unit  angular  velocity  El 

Vertical  force  per  unit  angular  displacement  Gl 
Pitching  moment  per  unit  angular  acceleration  A2 
Pitching  moment  per  unit  angular  velocity  B2 

Pitching  moment  per  unit  angular  displacement  C2 
Pitching  moment  per  unit  vertical  acceleration  D2 
Pitching  moment  per  unit  vertical  velocity  E2 


Pitching  moment  per  unit  vertical 
displacement 

Al/(0.  5.  RHO.  I3) 
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Gl/(0.5.  RHO.U2  I2) 
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4, 


B2/(0.5.RHO.U-  1  ) 
C2/(0.5.RHO.U2    i3) 
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G2/(0.5.  RHO.U2  l2) 
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C22 

D22 
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I.     INTRODUCTION 

The  unstable  motions  of  planing  craft  have  been  under  study  for 
many  years  and  have  been  the  subject  of  much  literature  (see  biblio- 
graphy). The  ability  to  be  able  to  predict  the  stability  characteristics 
of  a  particular  hull  In  the  early  stages  of  design  is  of  importance  to 
naval  architects.  A  knowledge  of  the  effects  of  variables  such  as  beam, 
deadrise  angle,   etc.  on  stability  would  permit  intelligent  corrective 
action  to  be  taken  to  increase  the  dynamic  stability  of  existing  craft. 

The  problem  of  planing  craft  stability  involves  many  variables 
and  empi  ricai  relations  between  some  of  the  design  variables  have 
been  developed  to  predict  dynamic  stability. 

Two  formulae  recently  developed  emperically  from  expfrimental 
data,  (2)  and  (12),   result  in  conflicting  design  requirements  to  increase 
stability  (see  Appendix  C).  It  is  therefore  desireable  to  develope  a 
theoretical  approach  to  the  problem  so  that  the  effects  of  design  var- 
iables can  be  determined  independently  of  experimental  data. 

Perring  (10)  attempted  a  theoretical  approach.  His  lack  of 
success  can  be  attributed  to  a  number  of  causes.   The  foremost  of 
these  being  lack  of  sufficient  experimental  and  theoretical  information 
to  predict  the  stability  derivatives  accurately  and  the  ommission  of 
Important  terms. 
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n.     THEORY 

A  planing  hull,  as  a  rigid  body,  has  six  degrees  of  freedom.  This 
study  treats  the  boat  as  a  two  degree  of  freedom  system  by  investigat- 
ing what  is  considered  to  be  the  most  important  motions,  heave  in  the 
Z -direction  and  pitch  about  the  Y-axis. 

The  equations  of  motion  are  nonlinear.  To  facilitate  the  solu- 
tion of  these  equations  it  is  necessary  to  linearize  them. 

Linearized  equations  of  motion  for  ships  have  been  developed 
by  Abkowitz  (1),   Korvin-Korvosky  (7)  and  others.  Those  of  Abkowitz 
are  most  complete.  If  the  coefficients,  i.e.  stability  derivatives,  are 
substituted  into  these  equations  and  then  the  result  transformed  into 
the  frequency  domain,  it  should  be  possible  to  evaluate  the  stability 
by  using  the  Routh  criterion  (5). 

The  method  used  to  predict  stability  or  lack  thereof  proceeds 
as  follows: 

1.  The  stability  derivatives  are  determined,   see  Appendix  A. 

2.  The  stability  derivatives  are  substituted  into  the  linearized  equations 
for  ship  motion.  (1) 

3.  The  resulting  equations  are  transformed  by  substitutions  of  the  form 
z  B  Zmax  e      and   8  ■    ©max  e    . 

4.  An  equation  in  S  is  obtained. 

5.  The  Routh  discriminant  is  evaluated  for  the  fourth  order  equation 
in  S,  see  Appendix  A. 
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III.     DESCRIPTION  OF  WORK  ACCOMPLISHED 

The  firststep  toward  a  solution  was  to  determine  the  stability 
derivatives  and  combine  them  to  form  the  coefficients  of  the  linearized 
equations  of  motion,  see  Appendix  A.  The  coefficients  were  then  non- 
dimensionalized  using  beam  as  the  non-dimenstonalizing  length  (12). 

A  computer  program  was  written  to  solve  for  the  Routh  discrim- 
inants, see  Appendix  H,  using  as  input  data  the  results  from  a  series 
of  tests  run  at  the  David  Taylor  Model  Basin  (2),  see  Table  3.  The 
resulting  discriminants  were  then  plotted  against  speed  showing  a 
consistency  in  the  directions  of  the  paths,  however  there  was  no  ob- 
vious difference  between  stable  and  unstable    boats,  see  Figure  9. 

At  this  point,  an  attempt  was  made  to  determine  the  roots  of  the 
fourth  order  stability  equation  to  examine  their  loci  using  a  computer 
program  from  the  MIT  "SHARE"  library.  (SHARE  No.   1514  RTSCH). 
RTSCH  proved  to  be  unsatisfactory.  The  answers  obtained  from  this 
program  are  seriously  in  error  even  for  the  simplest  of  input  equa- 
tions . 

It  was  then  decided  to  vary  in  turn  what  seemed  to  be  the  most 
important  variables:  DIFB1,   DIFC1,   B2,   Dl,   D2,   El,  Gl,   and  G2. 
This  was  done  to  determine  the  effect  of  changes  in  their  magnitude 
on  the  Routh  discriminant.  From  Figure  10  it  can  be  seen  that  vary- 
ing G2  roughly  grouped  the  stable  and  unstable  boats  with  the  unstable 
group  centered  about  G2  equal  to  0.  38  x  G2  at  the  point  of  zero  Routh 
discriminant.  The  program  was  then  run  with  G2  equal  to  0.  38  x  G2 
so  that  the  loci  of  the  discriminants  could  be  examined.  The  results 
are  shown  in  Figure  11.  Based  on  these  results  it  was  decided  to  in- 

1.  Unstable  boats,  as  refered  to  in  this  paper,   are  those  which  porpoised 
at  a  Fv    less  than  6.  0.  Stable  boats  are  those  which  had  not  porpoised 
before  maximum  test  speed  (2)  was  attained  (Fv   *  6.  0). 
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vestigate  G2  further.  Since  the  term  G2  is  made  up  of  two  parts,  force 

5(arm)                     S(force)     ,*  .      .  .  .. 

x   __ +  arm  x    __ it  was  decided  to  vary  these  two  parts 

independently  to  see  if  better  correlation  could  be  achieved  in  either 
of  the  two  groups.     Correlation  was  not  improved,   see  Figures  12 
through  17.  However  it  was  observed  that  the  discriminants  for  the 
stable  boats  changed  very  little  with  changes  in  G2,   or  its  parts.  On 
the  other  hand  the  discriminants  of  the  unstable  boats  changed  a  great 
deal  with  changes  in  G2.  This  lead  to  the  conclusion  that  there  must  be 
a  term  in  the  coefficients  of  the  characteristic  equation  which  over- 
powered G2  when  the  boat  was  stable  but  was  of  the  same  magnitude 
as  G2  when  the  boat  was  unstable. 

Based  on  information  obtained  thus  far,   it  was  decided  to  investigate 
each  term  of  the  Routh  discriminant  to  see  which  of  the  coefficients  were 
controlling  for  the  stable  and  unstable  boats.   Values  of  the  coefficients 
obtained  from  program  1  in  Appendix  H  were  inserted  into  each  term 
manually  and  inspection  of  the  results  was  unfruitful.  No  obvious  dif- 
ference could  be  detected  between  stable  and  unstable  boats.  It  was 
concluded  that  the  interaction  was  much  more  subtle. 

A  closer  inspection  of  each  term  indicated  that  the  coefficients  Z. 

and  M  •  ,  Dl  and  D2,  may  interact  with  important  effects  and  this  became 
w 

the  final  step  in  the  investigation  of  the  coefficients.  The  results  are 
shown  in  Figures  18  through  20.  At  this  point  the  investigation  of  the 
coefficients  of  the  equations  of  motion  was  terminated  because  of  time 


1.  At  this  point  it  was  necessary  to  reduce  the  number  of  plots  to 
three  stable  and  three  unstable  boats  in  order  both  to  simplify  the  plots 
and  to  make  more  efficient  use  of  computer  time.  The  unstable  boats 
were  selected  by  choosing  two  which  had  axis  intercepts  fairly  close 
together  and  third  whose  intercept  was  remote  from  these  (models 
4665-3,   4666-13,   and  4668-9  in  Figure  10). 
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limitations . 

Concurrently  with  the  above  work,   a  program  was  written  to  solve 
for  all  the  hydrodynamic  performance  characteristics  of  the  planing 
hull:  planing  angle,   wetted  surface*  resistance,   power  requirements, 
and  stability.   This  program  uses  design  parameters  as  inputs  and 
provides  an  easily  readable  output.   Development  of  the  criterion  for 
the  equilibrium  planing  condition  is  shown  in  Appendix  E  and  details 
of  the  program  are  contained  in  Appendix  H.  In  order  to  facilitate  the 
writing  of  this  program,   it  was  necessary  to  determine  an  expression 
for  mean  bottom  velocity  based  on  an  emperical  function  of  deadrise 
angle,  Appendix  D. 

The  coefficients  of  the  equations  of  motion  and  the  Routh  discriminant 
based  on  the  computed  planing  conditions  were  compared  with  those 
based  on  experimental  data.  The  result  of  this  comparison  is  shown 
in  Table  2. 

The  program  was  then  run,  for  model  4668-9  with  G2  equal  to  0.  38 
x  G2,  varying  BETAI,   EPSILI,  VCG,  BEAM,   CG,   and  YI  in  turn.  The 
results  were  then  plotted,   Figure  21, so  that  a  comparison  of  the  re- 
lative effects  of  the  variables  could  be  made. 
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IV.     DISCUSSION  OF  RESULTS 

The  plots  of  the  discriminant  versus  speed  obtained  from  the  first 
calculations,  Figure  9,   are  disappointing.  According  to  these  results 
most  models  were  stable  throughout  the  entire  range  of  speeds  in- 
vestigated,  contrary  to  the  experimental  results.  The  lack  of  agree- 
ment could  be  caused  by  one,  or  both,  of  the  following: 

(a)  incorrect  formulation  of  one  or  more  of  the  stability  derivatives 

(b)  neglecting  the  cross  coupling  effects  of  longitudinal  motion. 

Perring  (10)  indicated  that  inclusion  of  the  longitudinal  motion  cross 
coupling  effects  had  negligible  effect  on  the  outcome  of  his  solution  to 
the  problem.  It  is  possible  that  in  the  present,   more  refined  solution, 
the  magnitude  of  this  cross  coupling  may  become  relevant. 

The  investigation  of  the  effects  of  varying  the  magnitudes  of  the 
stability  derivatives  indicates  that  a  solution  to  the  problem  may  lie 
in  this  area.  Although  variations  in  Z    ,   Z  ,   and  (Z    *  u0   Z    ),   DIFB1, 

Mr  4~t  XJ  lnr 

DIFC1,   and  Gl,  failed  to  yield  any  evidence  of  consistent  influence, 
Figure  10  shows  that  variation  of  G2  produced  a  fairly  consistent  dif- 
ference between  stable  and  unstable  boats.  It  is  true  that  there  is  con- 
siderable scatter  of  the  axis  intercepts  within  each  group,  but  there  is 
an  undeniable  consistency  in  the  grouping.  It  is  also  of  interest  to  note 
there  is  much  less  scatter  in  the  stable  group  than  in  the  unstable  group. 
The  grouping  of  the  stable  boats  cannot  be  attributed  to  their  equal  Froude 
number.  An  examination  of  the  unstable  boat  grouping,  Figure  10,  in- 
dicates that  models  4665-3  and  4668-9  intercept  the  axis  at  the  same 
point  with  Fv  of  3.24  and  5.  03  respectively  whereas  model  4666-17, 
which  has  a  Fv  of  5.  01  (essentially  equal  to  that  of  model  4668-9),  has 
an  intersection  remote  from  the  preceeding  two. 

The  investigation  of  the  loci  of  the  discriminants  with  G2  equal  to 
0.  38  x  G2,  Figure  11,   show  that  the  original  curves,   shown  in  Figure  9, 
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now  bend  down  toward  negative  values  of  discriminants  as  speed  in- 
creases. Some  of  the  unstable  models  have  loci  consistent  with  ex- 
perimental data,  i.e.  the  Houth  discriminant  heads  toward  negative 
values  as  the  porpoising  speed  is  approached  however,  the  boats  which 
were  stable  throughout  their  test  range  have  loci  lying  completely  below 
the  axis  {negative  discriminants  indicate  instability).  This  indicates 
that  simply  multiplying  G2  by  a  single  factor  does  not  produce  reliable 
Houth  discriminants. 

The  investigation  of  the  effects  of  varying  the  two  components  of  G2 
separately,  Figures  12  through  17,  although  not  producing  better  corre- 
lation within  the  stable  or  unstable  groups  does  point  out  the  small  ef- 
fect that  these  variations  have  on  the  intercepts  of  the  stable  group 
compared  to  their  effect  on  the  unstable  boats.  This  information,  as 
it  stands,  indicates  that  other  terms  in  the  equations  of  motion  are 
more  powerful  at  stable  speeds  but  that,  at  the  porpoising  speed,  G2 
is  a  powerful  term. 

The  results  of  the  investigation  of  the  simultaneous  variations  of  Dl 
and  D2,  see  Figures  18  through  20,  point  the  way  to  what  may  be  a 
valuable  area  for  further  study.  The  first  useful  bit  of  information  ob- 
tained is  the  fact  that  Z.,   Dl,  has  very  little  effect  on  the  magnitude 
of  the  Houth  discriminant  and  that  M  . ,  D2,  has  a  large  effect.  The 
most  important  result  of  this  investigation  is  the  fact  that  the  axis  in- 
tercepts of  models  4666-13  and  4668-9  have  been  reversed  in  their 
relative  positions  from  what  they  were  when  the  coefficient  G2  was 
varied.  This  means  that  a  simultaneous  variation  of  D2  and  G2  may 
cause  these  two  extreme  boats  to  cross  the  axis  at  the  same  point  and 
thereby  correlate  the  unstable  group. 

Correlating  the  results  in  this  manner  does  not  really  solve  the 
problem.  The  stability  indicator,   Routh  discriminant,  is  not  reliable, 
as  the  discussion  of  Figure  11,  Appendix  F,has  shown.  Further  work  is 
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required  to  produce  better  agreement  between  the  intercepts  of  the 
discriminants  and  experimental  data.  An  experimental  investigation 
of  the  individual  stability  derivatives  for  comparison  with  the  theoreti- 
cally developed  derivatives  would  be  helpful  in  locating  the  terms  of  the 
equations  of  motion  which  need  to  be  reevaluated. 

The  program,   which  solves  for  the  hydrodynamic  performance 
characteristics  of  planing  boats,  yields  information  of  importance  to 
design. 

Starting  with  attitude,   wetted  keel  length,   wetted  chine  length,   and 
drag,  it  can  be  seen,  Table  1,  that  there  is  good  agreement  between 
theory  and  experiment  for  boat  attitude  and  drag.  The  theoretical  values 
of  wetted  keel  length  and  wetted  chine  length  are  larger  than  the  ex- 
perimental values. 

A  plot  of  WKEEL  -  WCHINE  vs  TRIM  comparing  theory  (12)  and 
values  calculated  from  the  experimental  results  of  (2),  Figure  22,   in- 
dicates that  the  mean  line  of  data  points  lies  above  the  theoretical  line 
for  this  group  of  boats,   Table  3.  The  largest  errors  occur  at  the  largest 
angles  of  attack. 

The  expression  developed  for  mean  bottom  velocity,   Appendix  D, 
yields  results  which  compare  very  accurately  with  graphs  shown  in 
Figure  7. 

The  comparison  of  derivatives  calculated  directly  from  the  program 
shows  good  agreement  with  the  exception  of  C2  and  E2,   see  Table  2. 
This  error  was  most  likely  caused  by  the  difference  in  actual  and 
calculated  wetted  length. 

The  results  of  the  variation  of  BETAI,  EPSILI,   BEAM,   YI  and  CG, 
shown  on  Figure  21,  can  not  be  conclusive  because  of  inconsistencies 
which  have  been  found  in  the  discriminant.  However  Figure  21  does  show 
that  the  beam,   deadrise  angle  and  longitudinal  position  of  the  center  of 


■ 
■ 

I 

.  "  ' 

©nil  r 

o  /  srit 

;:ia 


-  17  - 


gravity  have  the  largest  effect  on  stability.  The  implication  of  Figure  21 
that  a  decrease  in  beam  increases  stability  agrees  with  (12).  The  in- 
ference that  an  increase  in  deadrise  angle  results  in  a  decrease  in 
stability  is,   at  first  glance,   distressing.  However  an  inspection  of 
Figure  16  of  reference  (12)  indicates  that  this  may  well  be  the  case 
for  the  boat  examined.  An  increase  in  deadrise  angle  results  in  an  in- 
crease in  trim  and  a  decrease  in  the  lift  coefficient.  Both  of  the  latter 
effects  are  destabilizing.  If  the  destabilizing  influence  caused  by  the 
change  of  trim  and  of  the  lift  coefficient  is  greater  than  the  stabilizing 
effect  of  the  increase  in  deadrise  angle,  the  boat  will  be  destabilized. 

The  results  Indicate  that  moving  the  center  of  gravity  forward  in- 
creases the  range  of  stability.  This  forward  movement  results  in  a 
decrease  in  trim  and  it  is  known  that  a  decrease  in  trim  results  in  an 
increase  in  the  range  of  stability. 
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V.     CONCLUSIONS 

Unfortunately  this  thesis  does  not  go  far  enough  to  solve  the  problem 
of  predicting  the  stability  characteristics  of  a  planing  craft  in  the  de- 
sign stage.  The  relative  magnitudes  of  the  force  and  moment  derivatives 
which  make  up  the  coefficients  of  the  equations  of  motion  give  rise  to 
inconsistencies  in  predicting  stability  which  have  not  been  resolved. 

This  study  does  indicate  however  that  variations  in  the  following 
quantities  have  a  minor  effect  on  the  magnitude  of  the  stability  in- 
dicator (Routh  discriminant): 

(a)  change  in  lift  coefficient  x  area  with  respect  to  vertical  velocity 
(DIFB1) 

(b)  change  in  lift  coefficient  x  area  with  respect  to  vertical  position 
(DIFC1) 

(c)  change  in  vertical  force  with  respect  to  the  angular  acceleration 
about  the  pitch  axis  (Z«  or  Dl) 

q 

(d)  change  in  vertical  force  with  respect  to  the  angular  velocity  about 
the  pitch  axis  (Z    +  u0.  VERM  or  El) 

(e)  change  in  vertical  force  with  respect  to  trim  (Z    +  u0*  Z     or  Gl) 

tf  w 

(f)  change  in  pitching  moment  with  respect  to  velocity  in  the  heave 

direction  (M    +  uD'  M  »  or  B2) 
q  w 

and  that  the  following  have  a  major  effect: 

(a)  change  in  pitching  moment  with  respect  to  vertical  acceleration 

(M  .  or  D2) 
w 

(b)  change  in  pitching  moment  with  respect  to  vertical  position 

(M    or  G2). 

z 

The  computer  program,   developed  as  part  of  this  thesis,   which  solves 
for  the  other  hydrodynamic  performance  characteristics  of  planing  craft 
is  able  to  reproduce  experimental  results  with  minor  limitations.  The 
expressions  used  in  computing  wetted  length  of  chine  and  keel  do  not 
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yield  accurate  results.   This  causes  an  error  in  wetted  length  and 
aspect  ratio  which  is  further  reflected  in  some  of  the  coefficients  of 
the  motion  equations. 

The  expression  developed  for  computing  the  mean  bottom  velocity 
yields  consistently  good  results. 

The  results  of  the  study  of  the  variations  of  design  parameters 
generated  by  means  of  computer  Program  2,  Figure  21,  show  that  the 
range  of  stability  may  be  increased  by  decreasing  deadrise  angle,  de- 
creasing beam  and  moving  the  center  of  gravity  forward.  Changes  in 
the  vertical  center  of  gravity,  moment  of  inertia  and  shaft  angle  have 
minor  effects  on  the  range  of  stability. 

In  spite  of  the  inconclusive  results  of  the  stability  investigation,  the  re 
is  substantial  indication  that  the  stability  problem  can  be  solved.  The 
solution  of  this  problem  will  utilize  results  like  those  shown  in  Figure 
21  in  conjunction  with  Program  2.  This  should  provide  an  extremely 
useful  design  tool  for  optimizing  planing  craft  design.  For  example; 
assume  that  it  is  desired  to  design  a  planing  hull  to  operate  at  40  knots. 
The  design  procedure  would  proceed  along  the  following  lines: 

(a) Run  computer  Program  2  for  a  number  of  combinations  of  beam, 
deadrise  angle,  and  the  longitudinal  position  of  the  center  of  gravity 
obtaining  drag  information  for  ail  combinations  which  yield  designs 
stable  to  40  knots. 

(b)  From  the  data  thus  obtained,  develop  a  family  of  curves  for  each 
deadrise  angle  by  plotting  drag  versus  beam  for  several  locations 
of  center  of  gravity. 

(c)  Choose  the  design  for  minimum  drag  for  a  40  knot  planing  hull  from 
the  curves. 


- 


u-r  tjB  staler 

. 

•n^h  w  Bcroite 

■ 


-  20  - 


VI.     RECOMMENDATIONS 

1.  Repeat  the  study  described  herein  using  three  degrees  of  freedom: 
pitch,  heave,  and  surge. 

2.  Develop  a  more  accurate  method  of  predicting  wetted  length  of 
keel  and  wetted  length  for  a  planing  surface  with  deadrtse  angle  based 
on  experimental  data. 

3.  Make  an  experimental  investigation  of  the  force  and  moment 
derivatives  for  comparison  with  those  developed  theoretically. 

4.  An  examination  of  the  loci  of  the  roots  of  the  characteristic  equa- 
tion in  the  S-piane  would  produce  valuable  results  once  the  inconsistencies 
in  predicting  stability  are  ironed  out.  An  investigation  of  this  sort  would 
show  the  effect  that  variations  in  design  parameters  have  on  how  the 
roots  approach  and  cross  the  imaginary  axis.  This  requires  a  more 
accurate  root  extraction  computer  program  than  was  available  to  the 
authors . 
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APPENDIX    A 

Equations  representing  force  and  moment  equilibrium  for  a  ship  can 
be  expressed  in  the  form:  (1) 

A1Z    +    BIZ    +    C1Z    +    D19    +    E19    +  G19    =    Felwt    forces 

A29    +    B20    +    C20    +    D2Z  +    E2Z  +  G2Z  =  Melwt moments 

For  the  case  under  consideration  here,   smooth  water,  there  are  no 

iwt  iwt 

extsrnal  excitations.  Therefore  Fe       and  Me       are  both  equal  to  zero. 

st  st 

By  substituting  the  transforms  Ze      and  9e      into  the  force  and  mo- 
ment equations  we  can  obtain  a  characteristic  equation  of  the  form: 
AA.  S*+  BB.  S3  +  CC.  S2  +  DD.  S  +  EE  =  0 
Where: 
AA  =  1. 

A22.  Bll  +  All.  B22  -  D22.  Ell  -  E22.  Dll 


BB  ■ 

CC  = 

DD  • 

EE  = 


AU.A22  -  D11.D22 

A22.C11  +  B22.B11  +  A11.C22  -  D22.G11  -  E22.E11  -  G22.D11 

All    A22  -  Dll  D22 

B22.C11  +  B11.C22  -  E22.G11  -  G22.E11 
A11.A22  -  D11.D22 

C22.C11  -  G22.G11 
A11.A22  -  D11.D2~2 


If  the  Routh  criterion  is  applied  to  the  characteristic  equation  it  should 
be  possible  to  evaluate  the  stability  of  the  boat  (5).  The  criterion  indi- 
cates that  a  boat  will  be  stable  and  nonosciilatory  in  the  steady  state  if: 

BB.  CC.  DD  -  AA.  DD2  -  BB2.  EE  >  0 

Negative  Routh  discriminants  are  not  meaningful  for  the  case  under 
study.  Negative  roots  denote  instability,  instability  implies  motion,  and 
motion  in  this  case  implies  nonlinearity.  Since  the  method  is  based  on 
linearized  equations,  the  last  meaningful  Routh  discriminant  is  zero. 
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THE  DERIVATION  OF  COEFFICIENTS  FOR  FORCE  AND  MOMENT 
EQUATIONS: 


Al  and  A2  will  have  to  have  terms  containing  added  mass  and  added 
inertia  respectively.       This  is  necessary  because  the  added  mass  and 
added  inertia  may  be  of  the  same  order  of  magnitude  as  that  of  the 
boat  itself. 
Al  =  W  +  VERM 


Bl 


Bl  =  Z 
w 

Z         —  =  —■    *Z 
w       S  w         c  w       S  9 

a  9  ,   J^ 

*w  *    V 

h  Z  -     !      RHO    n*      *CLA 

Z     *  0.5    RHO    u0  DIFB1 
w 

Bl  *  Z 
w 

where: 

CLA  =  C-    A   (lift  coefficient  x  area) 

CLA 
DIFB1  -    -r-g —    (DIFB1  is  the  symbol  used  in  the  computer  program) 


CI 

CI  ■ 


I 


Z      .    _*    .    OM!tL=    1    rho.u!.     ?<CLA) 

z         d  z  ^z  2  °  &  • 


Z     *  0.5    RHO    u0a  DIFC1 
z 

CI    »   Z 
Z 


See  Appendix  B. 
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where: 


DIFC1  =  computer  program  symbol  for 


3CLA 
dz 


Dl 


Dl  =  Z. 

q 

q        Sq 


r  = 


VCG 


sin  a 
Limits  on  a 


e,%1 


F/GURE  I 


■l.CG-WETL,        7T  -1  .  CG   . 

COt       (— VCG >-2"*Cot       (VCG) 

VCG.    q 

v   =  r.   q  =    — : — 

sin  a 

VCG.    a.    cos(  0  +  a) 

v     *     - — : 

z  sin  a 

d  Z  =  d  (m.   w)  -  4"  •  RHO.  BEAM2,   d  £  (  ^5  ,   cos(  0  +  a  )  )  a* 

e>  sin  a 

I  =  VCG  .   cot  a 


See  reference  (9),    page  420,    Fig.    62A. 
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d  i  =    -VCG.    csc2ada 


\  =  CG 

Z_  =    J  .   rho  .   BEAM3.   VCG3|         cos  (9  + a)  da 

J  sin     a 

\  =  CG  -  WETL 


Z.   ■    ? 

q      q 


7T 


Z.  =    7T-.   RHO  .   BEAM3.   VCG  .   WETL; 
q      2 


LVCG(CG--2— )"Sine 


Dl  =  Z. 

q 


El 


0,1,1 


FIGURE  Z 


El  =  Z     +  u0.    Z . 
q  w 

sz      dw     se     sz        a  w    _, 

Z     =  - —   * .     - —  .   -r-r—    * .    Bl 


q       Sq         c  q'      d  w  *    B  9 


a  q 


We  now  wish  to  express  q  as  an  effective  velocity,    w. 
VCG 


r  ■ 


sin  a 


v  *  r  .   q  .   cos  (a   +  9)  ■  VCG  .   q  (cos  9  .   cot  a    -  sin  9) 


d  (v.  a)  *    increment  of  velocity  x  area 

=    VCG  .   q  .   (cos  9  .   cot  a    -  sin  9)  d  £  .   BEAM 
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V.A  «    -VCG2.   BEAM  .   q.fcos  9  J  cot  a    .    csc3ad  a    -  sin  0jcsc2a  da] 

where:  A  *  wetted  area  s  WETL  .   BEAM 

V.A  =  -WETL  .   BEAM  .   q.Icos  0  (0.5  WETL  -  CG)  +  VCG  .   sin  el 

Z     »  Bl  [cos  0  (CG  -  0.  5  WETL)  -  VCG  .   sin  0~| 

El  «  Z  ■  +  u0  .  VERM 


i*  CG 

4=CG-WETL 


Gl 


Gl  »  Zn  +  u0.   Z 

0  W 


Lift  =  -r-.   RHO  .   u2  .   CLA  =  Z 


*Z         1_  a  B(CLA) 

dT        2    °   RH°  '   U°       3  0 

ze  *  T  e  RHO  •  u°  DIFB1  3  uo  B1 

Gl  *  Z    +  u0  Z      ■  u0    Bl  +  u0   Bl 
w  w 

=  2.     u0  .   Bl 


A2    =  YI  +  VERYI* 

where: 

YI  ■  pitching  moment  of  inertia 

VERYI  =  added  inertia. 


See  Appendix  B  for  development  of  VERYI. 
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B2 


e,<?,<? 


FIGURE  3 


B2  = 

M     +  u0  . 

q 

M. 
w 

o 

M     = 

q 

d  (force 
dq 

.    arm) 

=  arm 

9  (force) 
a  q 

+  force 

/"aq 

f  (car  -  r 

PT.1  rns 

ft  -  vr< 

"1        sin   ft 

i    3  Z 

5q 


=  L(CG  -  CPL)  cos  9  -  VCG  .    sin  0  ]  .    El 


B2  =  M     +  u0  M  .    =  M     +  u0  .    D2 
q  w  q 


C2 


e,%1 


FIGURE* 
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r  =  (VCG2  -  CGS  ) 


.    -1 


WK  =  H/sin  9    ,      7    =  sin       (VCG/r) 


H  =  r.sin  (?  +  0)  -  (VCG  -  H        ) 

y  —  (j 

-  d(arm)  B  (force) 

M9  =  force   39 +  arm    —39- 

=  W    5  la5m    +  [  (CG  -  CPL)  cos  9  -  VCG  .  sin  9]  .Gl 

O  b 

Assuming:    CPL  =  Cx-WK 
B  (arm) 


d0 

MWK) 
d0 


d  (WK) 
=  -  CG .  sin  0  -  VCG  .cos  0  -  Cx.cos  0        *"        +  CX.WK  .  sin  0 


30 


=  -  r «  sin  ^  •  esc2  0  +  VCG  •  cot  0»csc  0 


=  VCG. esc  0  (cot  0    -  esc  0) 

\am)  =  CPL  [  sin  0  -  Y£§-  .cot  0  (cot  0  -  esc  0)  ]    -  CG  .sin  0  -  VCG  .  cos  0 

o  b  WK 

C2  =  M     +  u0    M      =  M.  +  u0  E2 
W  w  0 


D2 


FIGURE    5 


D2  = 

M. 
w 

M. 
w 

d  (force  .arm) 

d  w 


9  (force) 

arm  — - — ; +  force 

3   w 
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dM  =  f  .  RHO  .  BEAM2.  d£  (£.cos  0  -  VCG.sin  0).w 

Note:  Incremental  force  =  d  (added  mass  .  w) 

added  mass  term  taken  from  reference  (11),    page  420, 
Figure  62A,    for  an  increment  of  length. 

/CG 
Mw  =  T"  •rho*BEAm2J   (4- cos  9  -  VCG.sin  0)  d  4 

CG  -  WETL 
=   £.RHO  .  BEAM3. WETL  [cos  0  (CG  -  .  5  WETL)  -  VCG  .  sin  0  ] 

D2  =  M. 
w 


E2 


vcg.  sine 


FIGURE  6 


E2  =  M 
w 

„           SM          a  (force,  arm) 

w       a  w                  d  w 

P 

S(force)       e           d  (a 
=    arm  — r +  force  — -/ 

/m) 

a  w 


w 


=  ^(CG  -  CPL)  cos  0    -  VCG.sin  0  ] 


a  (lift) 
a  w 


M      =  [(CG  -  CPL)  cos  0  -  VCG.sin  0  ]   Bl 
w 


E2  =  M 


w 
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G2_ 

See  Figure  4. 
G2  =  M 


M 


z 

cM         ?j (force,  arm) 

SS  mi  2  '     ■  '" 

Z  Q  Z  nZ 

-            $  (arm)  d  (force) 

■  force   — r —  +    arm 


oz  Sz 

B(arm)  m       m  m  .     3  (CPL) 

*      -  COS  8 


Sz  B  z 

assuming:       CPL  «  CX.WK 
where  C,,  is  an  arbitrary  coefficient. 

WK  =     ,H  •% 
sin  6 

where   H  *  draft  at  transom. 

(WK)  1       5H  1 

a  z  sin  0    S  z  sin  8 

S(arm)  cos  9 


s        - 


d  z  sin  0 

CPL 
WK 


Cx  *    -  Cxcot  0 
cot  0 


PPT 

M    *    L(CG  -  CPL)  cos  0  -  VCG  sin  0  ]  Cx  -  W  ^~£   cot  0 
z  WK 

G2  =  M 

I 

Note:  Computer  program  determines     — by  incrementing  variables 

c  z 

in  the  equation   CPL  =  0.  75    -  r-^- 


5.21(-~)a  +    2.39 


developed  in  reference  (10),  page  16. 


=  SO 


M 


I 
M 

(J  «  (mi&)6 

- —  9  boo  -    -    — 

5iV  10      :anhntr«*t 

6  a 

C3W)  . 

Is  9  nla 

ft  AM  -»  -    -     - 

TC 

-,D  [fl  nh  -  G  boo  (J^D  -  OD)  j    ■     M 

s 

M  =  SO 

■ 

ev  gn  nfiigo-xq -is*uq. 

1  .0  =  J^D   noitaups  sift  nJ 


ee.i  iS.e 


.81  d^Aq  ,(CI)  ©oro-xs^i  nl  bsqols 
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APPENDIX    B 

Solving  for  Added  Mass: 

Because  of  a  lack  of  information  about  deadrise  surfaces  the  added 
mass  is  calculated  as  for  a  submerged  elliptic  cylinder.  The  result 
will  be  divided  in  half  and  the  minor  axes  of  the  ellipse  will  be  set  to 
zero.  This  is  assumed  to  be  a  suitable  approximation  of  the  added 
mass  of  a  deadrise  hull  at  the  water  surface. 

For  flow  past  an  elliptic  cylinder  (page  251,  (8): 


A_        B 


w  *    -3-    +    7s 


where:        w  ■  complex  potential 

It; 
r»e    \  a  unit  circle 

A  ■  u  (b  cos  a  +  i  a  sin  o  ) 

B  »i  o>(aa-ba) 
4 

For  the  case  in  question  where  only  the  vertical  oscillation  is  being 

considered, 

w  *  0 

so: 


\tr 

A 

a  — — 

s 

U  (b 

cos 

a  +  i 

a 

sin 

a) 

W 

S 

w 

=  Ar 

dw 

•  A 

T  * 

1 
4 

AA 
? 

dr 

Now:         T»    ■   ?  IC|       -^rd5"  T  *  kinetic  energy 

*    -    T   i  p  12  w  i  (the  residue  AA)  ] 
4 

■    -    ~   i  P  2ir  i  U8(b8cos8  a  +  a^in^  ) 
4 

U8 
For  this  case  a   *  90°,  b  *  0  and  T    =     g-  p  ?r  a8 

but  T  *  y-  M  Us 

then  2  VERM  *  p  *  a8  per  unit  width 

and        VERM  ■  tt   (WETL)8(BEAM) 


a 

:aa*M  bsbbA  i 

,B  <zdi  B9ob\iu9  saHb*9b  iuoda  aoiiaraiotal  i  b  to  ©8i 

;1IJ9  bsgierndua  4  lol  a«  b« 
of  }«•  m  fl9  grfj  lo  abXJB  ncnlm  srii  bna  Mb;  v 

b'-bi*-6  ad  ^jsmixoTiqqB  sfdsifcra  *  ed  oi  b 

.»o*>ii;a  isJaw  adJ  i&  .am 

!  5S  ^$aq)  i9bnll-^o  oilqilfa  nfi  taaq  wott  i 

H  A 

=  w 

Iaitn9toq  X9iqrnco  *  :9i9rfw 

9l5  ■  1 

(  j>  nla  a  i  +  »  a  U  ■  A 

(•d  -"«>'<•  r- 

i 

yiisd  el  nolralUDeo 

-  w 
U         A 


tr>Ma  Dll9fli^  -   T  ^  j    I     .     p 

olasT  9rfi)  MSJ    q]    i    -     . 

§ 

( x^nli^a  +  &  *aoD&d)"U  \  ^2  *,  i   ~    .    . 
*  *  1  ***  0  «  d   ,-09  »    d  98*..  fcJri* 

•U  M 
riJbiw  tino  i»q  "a  *  o,  ■  MH  »d» 

<MAaST)«(JTHW)   p  -  m  baM 
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The  Added  Inertia: 

The  problem  here  is  to  determine  the  added  inertia  for  a  body 
rotating  about  some  point  other  than  its  center. 

For  rotation  about  a  point  other  than  the  center  of  the  body  the 
complex  potential  is  given  by: 

w  -  i  w  Zn  Z 
In  this  case  Z0  *  Z  because  it  is  on  the  real  axis. 
Z  »  a  cos  Tj+  ib  sin  r, 

then 

A        B 
w  *  -=-  +  ==Ta-iu>Zc(a  cos  ?j  +  i  b  sin  fj ) 

w  *  AT+    BT  s  +  i  w  Z0  a  cos  r.  +  b  «  Z0  sin  rt 

~  *  A    +  2BT  -  i  w  Z0a  sin  ^|+bwZ0  cos  r,  fe, 

y *  e 

<*r_  -    Ui       i 

dTj  ir 

^-*  A  +  2BT  +  -§2   (i-  b  cos  v]  -  a  sin  r,) 

do  T  i 

.  ^.A  ^B  __ 

wjiv  ■  ___   +   _.  -  i  A  w  Z0  (a  cos  t,  +  i  b  sin  :  ) 


—       2BB 


+  2BA+  ~^r-  -  i  2  BTw  Z0  (a  cos  T;+  i  b  sin  r, ) 


1W8Z08  ,  t  .  v  /    1       u  ,  X 

__-2-  (a  cos  tj  +  i  b  sin  r( )  ( —  b  cos  ?,  -  a  sin  r, ) 

T  «    -  —  i  o/wdw 

-    -    i   i  pTT     ~+2BB  +  2BA+^-dr 

4         -4c)  r  7  * 


(a  cos  t]  +ib  sin  t,  )  (  -b  cos  t,   -  a  sin  r} )  dT 


Ac)    -^  J 

+     f     -lAuZ.(a  cos  >,+  i  b  sin  %  )  djf 
y(c) 

►    /       -i2Bru  Z0(aco8Tj+ib8in?i)d5'J 
J(c) 


■  Fl  +  F2  +  F3  +  F4 


sbbs  i 


: 


I 

m 


fb 


,T   M 


1  C 

,       r-  7  AA 

nit  d  ,  3)0S^aSj- 

i 


^b  K  +  * 


lb  (t  nlB  i  -   ;  Be  d  i  +  *) 


■*  c, 


lb  (    .  n)«  d  SuA 

L  ^  7.  esi-A- 


: 
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Fl  «  2  it  i  (AA  +  2  BB) 
F2,2,i^Z02ab!/ii5!^!l  tSs£l  dr 

taut  T-  s  it] 

dj~*  i  elr'  limits:  0  -     2  * 

o 


F2  =  2  ?r  iwaZc8  ab. 

F3  *  i  A  a)  Zc  I       (a  cos  r,+  ib  sin  r?)  (cos  r/  +  i  sin  ?))  d  q 

■  i  jr  A  w  Z0  (a  -  b) 

fit 

F4  *  -  i  2B  o>  Z0   j       (a  cos  fj  +  i  b  sin  ., )  (cos  2  /.+  1  sin  2  q )  d  r, 

*  -  i  2B  <*>  Z0(o) 

*  0 

T     •  *  4-  iifwdw 

»  -    j   i  p  (2  *  i  (AA  ♦  2  BB)  +  2  *r  i  waZ0a  a  b 
+  iriAuZ0(a-b) 

In  this  case  A  ■  A  *  0  because  there  is  no  Z  translation. 

Therefore: 

T     *  -   i  i  p  (2  jr  i  (2  BB)  +  2  n  i  w2Z0aa  b) 

*  |(2ws  (aa+  b8)  ♦  u)3  Zn8  a  b) 
for  the  plate,  b  *  0 

so    T  *  -jr  <oa  a 

1  4 

VERYI    =  -r-irpa      per  unit  width 

*  4"  *  P  (WETL)4  BEAM. 
8 


r.b 


:   I 


N 


I 

j  d  1  *) 


a) 

%  i 

i  ■  -. 

3 

0  =  A 

riT 

■ 
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APPENDIX    C. 

From  Figure  16  of  (12)  the  stability  of  a  planing  craft  can  be  in- 
creased by  increasing 

i.e.     C. 


V     2  L' 

But   CL  «  .0120  X  l'2  t1,  l        (page  10,  (12)  ) 

Since  increasing  C.    increases  stability,   increasing  X  increases  stability. 

1  m 
In  this  case   X  *  — =-*    .  This  means  that  a  decrease  in  b  will  result  in  an 

b 

increase  in  stability. 


From  Figure  20  of  (2)  the  stability  criterion  is 

CLb  1.80 

cp'  (  Fv  ) 

where: 

W 

Q  =    II 

Lb        i-pVsba 

and:  . 

V        ~ 

f      =  <-JT>2 

«V3 
c 

If     - *  -? decreases  the  boat 

^cjp_        kv3b    1 
b  2  cp 

is  stabilized.  Therefore  increasing  b  stabilizes  the  boat. 

A  comparison  of  the  two  methods  makes  it  difficult  to  decide  what 
effect  b  has  on  stability. 


aj»<St)H 


fit  iST  ilJ  98*0109b  B   iMtit  B: 


al  noli©  ;IMlda*a  edf  <S)  t< 

.1 

) 

taod  9ll*  a»B.6^T39b      II 

qo 
.taod  stt  aasilldaia  d  jniaa 

illlb  11  6  >wJ  ari:,  A 

'lldstc 
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APPENDIX    D, 

In  order  to  determine  the  resistance  of  a  planing  surface  it  is  nec- 
essary to  know  the  mean  velocity  over  the  bottom. 

In  Figure  12  of  (12)  a  method  is  provided  for  determining  this  graph- 
ically. Part  of  an  analytic  solution  is  provided  which,  if  completed, 
would  be  useful  in  the  computation  of  mean  velocity  in  a  computer 
program. 

The  following  is  a  resume  of  the  development  of  the  complete  equation. 
From  (12): 

Y».      L  0  0120^-1f~ 

V  X  2*  COS  T 

where: 

VM  «  mean  velocity  over  bottom  (fps) 

V     *  the  u0  ,  the  horizontal  velocity  of  the  origin  of  coordinates  (fps) 

t      *  trim  angle  of  planing  surface  (degrees) 

f((3)  =  an  undetermined  function. 

By  changing  1  to  radians  and  plotting 

LI 

.0120  r 
it  was  possible  to  find  an  f(p)  which  yielded  satisfactory  results. 
The  final  result  is: 


VM  I  0.  120  rlml  80.    -     50.  p 

V       "      J  XICOST  ■         ~^T- 

Wheref(9>-    80-    -B0-P 

COS      p 

t      *    trim  angle  in  radians 

f3      a    deadrise  angle  in  radians. 
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APPENDIX    E. 


EQUILIBRIUM  PLANING  CONDITIONS. 


FIGURE  8 


Note:  This  development  follows  reference  (12)  in  general.   However,   the 

final  result  contains  terms  neglected  by  SAVITSKY. 
Summation  of  forces  in  vertical  direction: 

(1)  W    =    N.cosG    +  T  .  sin  (0  +  €  )  -  Df.sin  0 
Summation  of  pitching  moments: 

(2)  N.c    +    Df.a    -    T.  f    =    0 
Summation  of  forces  along  keel: 

(3)  T.  cos  €  =  D     +    W.sin  0 


Combining  (1)  and  (3): 

W 
(4)    N    = 


D, 


COS    0. COS    € 


)S   c 


cose  -  sin  0.  sin  (0  +  e)J  +c^q^ 

[sin  0.  cos  e   -  sin  (0  +  € 


>] 
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Combining  (2)  and  (3): 


(5)    N.c    +  D..a 

I  COS   € 


W .  sin  9   +    D, 


=    0 


Combining  (4)  and  (5),  we  obtain  the  general  equilibrium  requirement: 
[  cos  c    -  sin  9    sin  (9  +  c)  j    c      -     f.  sin  9 


W 


cos  9  .  cos  e 


cos  c 


+  D. 


L  sin  9    cos  €  -  sin  (9  +  c)l  m  ,  .       t     1 

' — — c  +  a *  u 

cos  9      cos  €  COS  €J 


-  VE 


(«) 


:tne  oD 


9)  nia    G   . 


3    8C 


W 


9  n 


-  33  - 


APPENDIX     F 
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& *& 


+ 


* — * 

+ X 


SHAFT  ANGLE  /NCR£-AS£D  F^ATOM    7.  3°  ro  12.  5° 

VCG    ^DECREASED   FROM     0-3  FT      TO   0-2FT 

BEAM  DECREASED  FROM      /./3FT    TO     O&FT 

BEAM  INCREASED    FROM      /.IS FT  TO      /.-?  FT 

1 
^DEADR/SE"  DECREASED    F~ROM  J2-S     TO    S 


O-- 


fADR/SET  ///CREASED  FROM    /2-S    TO  20 


..!... 


OEADR/S. 
Vf  DECR 
YI   INCREASED     FROM      //.7 FT /.S  SECa  TO  20.0  FT ^S  SECe 


I  EASE/}    FROM    II.  7  F~r  LB  SEC*    TO  S.O  FTEB  SEC 


* 


„, 


E) tJ 

E S3  CC/HCR£-AS££>     FROM    2.35 FT    TO     S.SFT 


CS  DECREAS£D  FROM    2.S>S  F~T     TO   2.  S  F~T 


~ ' 





<?* 


2£ 


23 


SO 


32  3+  3*>  38  -+0 

VELOC/TY  (  ft /sec) 


+2 


■4~* 


+6 


I : I I I ! 


_. 
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Experimental  Results  from  (2) 


As  computed 
from  program  2 

Experimental 
Results 

Model  No.                     4668                          AT                 32.7     FPS 

Trim  (degrees) 

3.  68 

3.  80 

W  Keel  (ft.  ) 

5.09 

4.40 

W  Chine  (ft.  ) 

3.  78 

3.  10 

Drag  (lbs.  ) 

25.  73 

25.09 

Table  1.     Comparison  of  calculated  and  experimental 
hydrodynamic  characteristics. 
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Data 

For  Model  4668  Test  No.    9  from  (2)  for  the  19.  36 

Knot  Speed. 


Derivative 

As  computed  from 
experimental  planing 
conditions 

As  computed  from 
calculated  planing 
conditions 

All 

10.  16 

13.  16 

Bll 

3.  13 

3.  30 

Cll 

2.  43 

2.  38 

Dll 

6.  68 

6.  32 

Ell 

27.  44 

36.45 

Gil 

20.  44 

21.  55 

A22 

7.  25 

9.  47 

B22 

22.  61 

20.  58 

C22 

5.  59 

-1.  06 

D22 

6.  68 

6.  32 

E22 

0.  89 

-0.  14 

G22 

1.  20 

1.  20 

RD 

-0. 071 

0.  209 

Table  2.     Comparison  of  coefficients  of  equations  of 

motion  as  computed  from  experimental  data 
and  data  generated  from  computer  program  2. 
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Table  3.     Model  Data  (2).     All  models  have  deadrise  =  12.  5" 


Model 
Number 

Run 

Weight  (lbs) 

Beam  (ft) 

Length  (ft) 

CG  (ft) 

F       (at  max. 
test  speed) 

4665 

1 

54.  50 

1.  654 

3.912 

1.  62 

5.98 

4665 

3 

129. 08 

1.  654 

3.  912 

1.  70 

3.24 

4665 

7 

80.  07 

1.  654 

3.  912 

1.  70 

6.05 

4665 

8 

80.  07 

1.  654 

3.  912 

1.  55 

3.  50 

4665 

9 

80.  07 

1.  654 

3.912 

1.  39 

2.  74 

4665 

10 

55.  77 

1.  654 

3.912 

1.  86 

5.96 

4665 

11 

54.  50 

1.  654 

3.912 

1.  70 

-    5.99 

4665 

12 

54.  50 

1.  654 

3.912 

1.  55 

5.98 

4665 

13 

54.50 

1.  654 

3.  912 

1.  39 

3.  23 

4666 

9 

146. 20 

1.  623 

5.987 

2.  17 

3.  75 

4666 

13 

101. 80 

1.  623 

5.  987 

2.  17 

4.  53 

4666 

17 

76.  10 

1.  623 

5.  987 

2.  17 

5.  01 

4667-1 

9 

221. 10 

1.  600 

8.  00 

2.  95 

4.  02 

. -.— — — 

4668 

9 

141. 80 

1.  190 

8.  00 

2.  95 

5.  03 

4669 

1  6 

51.  40 

0.935 

8.  00 

3.27 

6.02 

_ -     ■   ■■  ■ 

Note:  Unstable  boats,    as  refered  to  in  this  paper,    are  those  which 
porpoised  at  F      less  than  6.0.    Stable  boats  are  those  which 
had  not  porpoised  before  maximum  test  speed  of  reference  (2) 
was  attained  (F     «  6.  0). 
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PROGRAM  1 


Program  1  can  be  used  to  solve  for  the  stability  derivatives,   co- 
efficients of  the  equations  of  motion,  the  coefficients  of  the  character- 
istic equation  and  the  Routh  discriminant  using  experimental  data  as 
input. 

Data  cards  are  punched  in  the  manner  indicated  by  READ  1  and 
1  FORMAT  where: 


W 

ALFAO 

CG 

C 

RT 

WK 

WC 

S 

TRIM 

YI 

RHO 

U 

BEAM 

BETAI 

VCG 

EPSIL 


weight  of  boat  in  (lbs) 

trim  angle  when  boat  is  at  rest  (degrees) 

longitudinal  position  of  center  of  gravity  forward  of 

transom  (ft) 

forward  speed  (knts) 

towing  force  (lbs) 

wetted  length  of  keel  (ft) 

wetted  length  of  chine  (ft) 

wetted  area  (fta) 

(planing  angle  -  a0)  (degrees)  is  change  of  trim  from 

the  at  rest  position 

moment  of  inertia  about  the  Y-axis.  Axis  taken  through 

center  of  gravity,  (lb  ft  sec3) 

density  of  water  (lb  seca/ft) 

arbitrary  non  dimensionalizing  velocity  (ft/ sec) 

beam  of  boat  (ft) 

deadrise  angle  (degrees) 

height  of  center  of  gravity  above  the  keel  (ft) 

shaft  angle  (degrees). 


o&imtio  ?<ci3  '      : 


btiM  1 


^riw  1  )1  | 


(■  b)  i**t  is  ex 

lo  b-iswiol 

(tl)  moanai* 

to  di%nsl  bettsw 
udo  to  rfigneJ  bettsw 

(8tl)  J891S  t 

«  li     (893  I§  3b  ;infl  gr 

risooiritns.vfii  BixA   ,8ixjB.y  9jf. 

(^\Sd8  di)  ir  Jan»b 

a\rt)> 

(tf )  *~  mudd 

(8»9«XJ£  a  »«ftU 

:/*i8  lo  i»tneo  lo  trisisri 

.(89r 


W 
OA" 

OO 

D 

I 

IY 

n 
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COMPUTER     PROGRAM    ONE 


C 

c 
c 
c 

c 
c 

c 
c 


c 
c 
c 


EQUATIONS  OF  MOTION  OF  PLANING  HULLS 

RD=ROUTH  DISCRIMINANT 

AA,BB,ETC=ROUTH  CRITERION  FACTORS 

All iBll »ETC=NONDIMENSIONAL  FORCE  AND  MOMENT  COEFFICIENTS 

Al tBl »ETC=FORCE  AND  MOMENT  COEFFICIENTS 

READ  INPt  T  DATA 

666  READ1 »W»Al  FAO , CG »C ♦ RT , WK ♦ WC »S * TR I M , Yl ♦  RHO, U ♦ BEAM ♦ BETA  I  , 
1VCG»EPSIL 
1  FORMAT(F7.2»F5.2,F4.2,F5.2,F5.2»F4.2,F4.2,F5.2,F4.2,F5,2. 
1F5.3*F3.0»F5.3»F4.1,F4.2>F5«2) 

V«C*1.689 

WETL= { WK+WC ) /2.0 

VERM=.125*RH0*3.1416*(WETL*^2 )*BEAM 

VERYI=.06  2  5*.062  5*RHO*3.1416*( WETL**4 )*BEAM 

CV=V/SORTF(32.2*BEAM) 

A=WETL/BEAM 

CPL=WETL*(0.7  5-1.0/( 5.21*1 ( CV/A ) **2 ) +2 . 3  9 ) ) 

TAU=(ALFAO+TRIM) / 5 7. 2  956 

BETA=BETAI/57.2956 

EPS=EPSlL/57.2956 

DIFFERENTIATION  OF  LIFT  COEFFICIENT  X  AREA  WITH  RESPECT 
TO  TAU 
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C 
C 


TAU1=TAU-0.001 

TAU2=TAU+n#001 

A=l./A 

CL VOL  1  =  1./ ( 2.*WETL* ( CV**2 )  )*(  (  (WC**2 )*SINF( 2.*TAU1 ) 
1/BEAM+1./3.*(2.*WC+WK)*SINF(BETA) /COSF(BETA) ) ) 
3  CLB1=0.5*CLV0L1 

CLSA1=(1.5708*A*TAU1*(C0SF(TAU1)**2)*( 1  .  0-S I NF ( BETA ) ) 
1/(1.0+A)+4#0*(SINF(TAU1)**2)*(COSF(TaU1)**3 )*COSF(BETA) / 
23.0+CLB1 )*S 

CLVOL2=l./(2.*WETL*(CV**2) )*( ( ( WC** 2 ) *S I NF ( 2 . *TAU2 ) 
1/BEAM+1./?.*(2.*WC+WK)*SINF(BETA) /COSF(BETA) ) ) 
5  CLB2=0.5*CLVOL2 

CLSA2=( 1.5708*A*TAU2*(COSF I TAU2 ) **2 ) *  I  1 . 0-S I NF ( BETA )  ) 
l/( 1.0+A)+4.0*(SINF(TAU2)**2 ) * ( COSF ( TAU2 ) **3 ) *COSF ( BETA ) / 
23.0+CLB2} *S 

CL*1.5708*TAU*A*COSF( TAU ) **2* ( 1 #-S I NF ( BETA ) ) / ( 1 ,+A ) +  4.* 
1SINF(TAU)**2*C0SF( T AU ) **3*COSF ( BETA ) / 3 . + ( CL VOLl+CL VOL 2 ) /4. 

TA0=TAU1 

TAB=TAU2 

CO=1.5  7  08*TAU*A*COSF(TAO)**?*( 1 . -S I NF ( 3  ETA )  )/ ( l.+A)+4.* 
lSINF(TAO)**2*COSF(TAO)**3*COSF(BET,\)/3.+(CLVOLl )/2. 

CB«1.5  708*TAU*A*COSF(TAB )**2*U.-SINF(BETA) )  /  Q.+A)+4.* 
LSI NF ( TAB )**2*COSF( TAB )**3*COSF< BETA) /3.+(CLVOL2)/2. 

DADTAU=10  00.*W/ (RhO*V**2 )*< l./CB-l./CO) 

DIFBl=(CLSA2-CLSAl)/.002+CL*DADTAU 

DIFFERENTIATION  OF  LIFT  COEFF  X  AREA  WITH  RESPECT  TO  Z 
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DELZ=0.001 

DELA=-DELZ*BEAM/SINF(TAU)/WETL**2 

DELS=BEAM*DELZ/(SINF( TAU ) *COSF ( BETA ) ) 

DELL=DELZ/SINF (TAU) 

CLV0L3=( ( (WC-DELL)**2)*SINF(TAU)/BLAM+(2.*WC+WK-3.*DELL) 
1*SINF(BETA)/(C0SF(BETA)*3.)  ) / I  2 . * ( WETL-DELL) * ( CV*#2 )  ) 
7  CLB3=0.5*CLVOL3 

CLSA3=(  l,fS708*(A-DELA)*TAU*(COSF(  TAU) **2 )*( 1 .-S INF ( BETA )) 
l/( l.+A-DELA)+4.*(SINF(TAU)*>2)*(C0SF( TAU )**3)*C0SF( BETA) /3. 
2+CLB3)*(S-DELS) 

CLV0L4=( ( (WC+DELL)**2)*SINF(TAU)/BEAM+<2«*WC+WK+3«*DELL) 
1 *SlNF( BETA) /(COS F(3ETA)*3. ) ) / ( 2 . * ( WETL+DELL) * ( CV**2 ) ) 
9  CLB4=0.5*CLVOL4 

CLSA4=( 1.5  708*(A+DELA)*TAU*(COSF( TAU) **2 ) * ( 1 . -S I NF ( BETA ) ) 
l/( l.+A+DELA)+4.*(SINF(TAU)**2)*(COSF( TAU ) ** 3 ) *COSF ( BETA ) /3. 
2+CLB4)*(S+DELS) 

DIFC1  =  (CLSA4-CLSA3) /♦002 

DIFFERENTIATION  OF  MOMENT  WITH  RESPECT  TO  Z 


CPLl  =  .7-5*  (WETL-DELL)  /(  5  .  2  1*  (  CV*BE  AM  )  **2  /(  WETL-DELL  )  **2  +  2  •  39  ) 

CPL2=. 7 5- (WETL+DELL) /( 5 . 21* ( CV*BEAM ) **2 / < WETL+DELL ) **2 +2 • 3  9 ) 

DCPLDZ  =  (CPL2-CPL1 ) /.002 

C1=0.5*RHO*( V**2 )*DIFC1 

DMDZ=(CG-CPL)*Cl/COSF(TAU)-W*DCPLDZ/COSF( TAU) 

WTCL=.5*RHO*S*V**2*CL 

Al=W/32.2+VERM 
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B1=0.5*RHO*V*DIFB1 

C1=0.5*RHO*(V**2)*DIFC1 

D1=VERM*VCG* (COSF( T AU ) * ( CG-. 5*WETL) /VCG-SINF( TAU) ) 

El=Bl*(C0SF(TAU)*(CG-WETL/2.)-VCG*SINF(TAU) ) 

G1=V*B1 

A2=YI+VERYI 

B2=El*( (CG-CPL)*COSF(TAU)-VCG*SINF(TAU) ) 

C2=G1*( <CG-CPL)*COSF( TAU ) -VCG*S I NF ( TaU ) ) +W* ( CPL* ( S I NF ( TAU ) -VCG* 
1C0SF(TAU)/SINF (TAU) /WK*(COSF( TAU ) /S I NF ( TAU ) -1 . /SI NF ( TAU ) ) )-CG* 
2SINF(TAU)-VCG*C0SF( TAU) ) 

D2=VERM*(COSF(TAU)*(CG-WETL/2. ) -VCG*S INF ( TAU ) ) 

E2=B1*( (CG-CPL )*COSF(TAU)-VCG*SINF(TAU) ) 

G2  =  DMDZ 

El  »  El  +  V*VERM 

G1=G1+V*B] 

B2  =  B2  +  V*D2 

C2=C2+V*E2 

QUANTITY  0.5  RHO  CANCELLED  OUT  OF  ALL  FOLLOWING 


A11=A1/(BEAM**3) 
B11=B1/(U*(BEAM**2) ) 
C11=C1/(BEAM*(U**2) ) 
D11=D1/(8FAM**4) 
E11=E1/(U*(BEAM**3) ) 
Gll=Gl/( (BEAM*U)**2 ) 
A22=A2/(BEAM**5  ) 
B22=B2/(t  *(3EAM**4) ) 
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C22=C2/( (BEAM**3)*(U**2) ) 
D22=D2/(BEAM**4) 
E22=E2/(U*(8EAM*»3 ) ) 
G22=G2/( (BEAM*U)**2 ) 
AA»1. 

BB=(A2  2*B11+A11*B22-D22*E11-E22*D11)/(A11*A22-D11*D22) 
CC= (A22*C11+B22*B11+A11*C22-D22*G11-E22*E11-G22*D11)/(A11 
1*A22-D11*D22) 
DD=(B22*C11+B11*C22-E2  2*G11-G2  2*E11)/<A11*A22-D11*D2  2) 
EE=(C22*C11-G22*G11 )/( Al 1*A22-D1 1*022 ) 
RD=BB*CC*DD-AA*(DD**2 )-(BB**2 )*EE 
PRINT39 

39  FORMAT  (  17X.3HA11»17X,3HBH»17X,3HC11»17X,3HD11»17X»3HE11. 

117X.3HG1]  ) 
PRINT  40fAlltBll»Cll»Dll»Ell.Gll 

40  F0RMAT( 6F20.5) 
PRINT  41 

41  FORMAT ( 15X,3HA22»15X,3HB22»15X,3HC12.15X,3HD22»15X,3HE22. 

115X,3HG22,15X»3H  RD) 
PRINT42»A22.B2  2»C2  2»D2  2»E2  2»G22»RD 

42  FORMAT( 7F3 8,5) 

PRINT43 

43  F0RMAT(14X,2HBB.14X,2HCC»14X,2HDD,14X.2HEE,13X,3HRD1.11X, 

15H0IFB1.11X,5HDIFC1 ) 
PRINT44»BB»CC»DDtEE»RDltDlFBlfDIFCl 

44  FORMAT(7F16.5) 
PRINTll»W.ALFAO»CG»V,RT,WK,WC*S»TRIM,YI  ♦  RHO  ♦  U  ,  BE AM ♦ BET A  I » 

1VCG.EPSIL 
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11    FORMAT(FT.2.F5.2,F4.2.F5.2,F5.2»F<r.2,F4.2.F5.2.F<f.2»F5.2» 
lF5.3.F3.0»F5.3,F<f.l,F4.2»F5.2//) 
GO    TO    666 
END 
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PROGRAM    2. 

This  program  solves  for  planing  conditions:  TRIM,  ASPECT  RATIO, 
WETTED  KEEL.,   WETTED  CHINE,    WETTED  AREA,   DRAG,   DRAFT  at 
the  transom,  MEAN  WETTED  LENGTH,  ESTIMATED  EHP  and  the 
STABILITY  INDICATOR  (Routh  discriminant:  positive  indicators  imply 
stability,  negative  indicators  imply  instability.   This  section  of  the 
program  does  not  yield  satisfactory  results  as  yet. ) 

A  listing  of  all  iterations  involved  in  solving  for  the  planing  conditions 
and  a  listing  of  the  coefficients  of  the  equations  of  motion  can  be  ob- 
tained as  indicated  by  comments  on  the  first  page  of  the  program  print 
out. 

This  program  uses  as  input,    1st  card: 
LIST  1,    LIST  2,   N  BOATS 
FORMAT  (3  I  3). 

LIST  1  and  LIST  2  are  as  defined  on  first  page  of  the  program  print 
out. 

N  BOATS  is  the  number  of  different  boats  to  be  run. 

2nd  card: 

BETAI,   EPSILI,   F,   VCG,   BEAM,   CG,   RHO,   YI.   W 
FORMAT    (4  F  5.  2,   5  F  10.  2) 

BETAI  ■  BETA    of  PROGRAM  1 

EPSILI  •  EPSIL  of  PROGRAM  1 

F  is  the  perpendicular  distance  from  shaft  center  line  to  CG  (ft). 
All  other  variables  are  as  defined  on  page  58. 

3rd  card: 

NUMBER,   IDENT 
FORMAT  (I  3,   5A4) 

NUMBER  *  number  of  speed  cards  which  are  to  follow 

IDENT  any  identifying  statement  or  symbol  not  to  exceed  20  spaces. 


»rf*  baA  ntU 

Xlqml  •;  3AT3 

9ti t  J  rlld*t8 

(  J;v  &ii  g  >io*iq 

anoiJibnoo  \  .,  &U  A 

o  ad  fl«o  rail  a  bna 

Jnliq  Bttlfetq  9^  *o  sgaq  til 

:b*xao  tel    .tuqni  sb  asau  aiJBi§oiq  »: 

.(-  hfUCI 

*iq  xtu  bsp  Isb  at 

:b*i*o  bnfi 

(S.OX  H  6   ,8  O^ 

I  MAHD 

alb  -u 

IMtW 
<*  MO? 

■ 
.a»o«qe  0£  beeoxs  o)  l  aI 
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4th  card: 

VKTS 
FORMAT  (F  10.  2) 

VKTS  is  velocity  of  boat  in  knots.   One  card  is  required  for 

each  speed.   Number  of  cards  equals  NUMBER  on  card  3, 


2nd   boat 


Figure  23.  Program  Assembly. 
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STABILITY  CHARACTERISTICS  FOR  PLANING  BOAT  SERIAL    6809 

EQUILIBRIUM  PLANING  CONDITIONS 

VELOCITY  (FPS)  =  25.  33 

TRIM  (DEG.)  ■  4.275 

ASPECT  RATIO  *  4.03 

WETTED  KEEL  (FT)  ■  5.  02  WETTED  CHINE  (FT)  «  4.  58 

WETTED  AREA  (FT**2)  *  5.  73  DRAG  (LBS)    p  21.  40 

DRAFT  (FT)  *  .37 

MEAN  WETTED  LENGTH  (FT)  =  4.  80 

ESTIMATED  EHP  •  .  99 

STABILITY  INDICATOR  -  .  13479E-01 


Figure  25.    Sample  of  Computer  Output  when  both  LIST  1  and 
LIST  2  equal  2. 
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C  COMPUTER  PROGRAM  TWO 

C      STABILITY  OF  PLANING  CRAFT 

DIMENSION  Tl(15)i  VALUE ( 1 5 )  ♦ I  DENT ( 5 ) 

COMMON  AA»  Allt  A22»  ASP.  BB.  Bill  B22,  BEAM.  BETA.  CC, 

1C11.  C22.  CLB.  CLO.  CG.  CV,  DD,  Dll.  D22.  DIFB1.  DIFC1. 

2DRAG.  EE.  Ell.  E22.  EPSIL,  F.  Gil.  G22.  LISTl.  LIST2.  NR. 

3RHO.  S.  TAU.  TRIM.  Tli  VALUE.  V.  VCG,  VKTS.  VM.  W,  WETL. 

4WCHINE.  WKEEL.  YI 
READ  1. LISTl.  LIST2.NBOATS 
C 

C      IF  LISTl*  1.  PRINT  OUT  ALL  COEFFICIENTS  AND  DERIVITIVES 
'C      ASSOCIATED  WITH  STABILITY  EQUATIONS 

C      IF  LIST  1=2.  PRINT  OUT  ONLY  STABILITY  INDICATOR 
C      IF  LIST  ;  =  1.  PRINT  OUT  ALL  ITERATIONS  INVOLVED  IN 
C      SOLVING  FOR  EQUILIBRIUM  PLANING  CONDITIONS 
C      IF  LIST  2=2.  PRINT  OUT  ONLY  FINAL  PLANING  CONDITIONS 
C 

1  FORMAT  (313) 

DO  7  NN=l.NBOATS 

READ2.BETAI  .EPSIL  I  ,F  .  VCG  .BEAM  »CG  .  RtiO,  Y  I  ,W 

2  FORMAT(4F5.2.5F10.2) 
BETA=BETAI/57.2956 
EPSIL  =  EPSILI/57.2956 

READ  3.NUMBER.  ( IDENT(  I  )  ,  1  =  1  .5) 

3  FORMAT( I3.5A4) 
PRINTlOlt  <  IDENTU  )  .1  =  1.5) 
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DO  7  M  ■  1  , NUMBER 

READ  4,  VKTS 

V  =  VKTS*  1,689 

4  FORMAT(F10.2 ) 
CALL  ANGLES 
CALL  C0EFF1 
EHP=V*DRAG/550. 
PRINT  116. EHP 

RD=BB»CC*DD-AA*DD**2-BB**2*EE 
IF(LIST1-1 )6»5,6 

5  PRINT102.AA. BB.CC.DD.EE 
PRINT103.A11 
PRINT104.B11 
PRINT105.C11 
PRINT106.D11 
PRINT107.E11 
PRINT108.G11 
PRINT109.A22 
PRINT110.B22 
PRINT111.C22 
PRINT112.D22 
PRINT113 « E22 
PRINT114.G22 

6  PRINT115.RD 

7  CONTINUE 

101  F0RMAT(43H1STABILITY  CHARACTERISTICS  FOR  PLANING  BOAT. 
19H  SERIAL   »5A4///) 

102  F0RMAT(46H  THE  COEFFICIENTS  OF  THE  FOURTH-ORDER  EQUATION. 
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153H  OF  MOTION  (AA*S**4  +  BB*S**3  +  CC*S**2  +  DD*S  +  EE  = 
2»33H  0.0)  ARE  AS  FOLLOWS  ( AA  THRU  EE ) /5F25.4// ) 

103  F0RMAT(46H  FORCE  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  VERT» 
119HICAL  ACCELERATION  =.F10.3) 

104  FORMAT(46H  FORCE  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  VERT, 
119HICAL  VELOCITY      =.F10.3) 

105  F0RMAT(46H  FORCE  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  VERT. 
119HICAL  POSITION      =.F10.3) 

106  F0RMAT(46H  FORCE  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  ANGU» 
119HLAR  ACCELERATION   =»F10.3) 

107  F0RMATC46H  FORCE  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  ANGU. 
119HLAR  VELOCITY       =.F10.3) 

108  F0RMAT(46H  FORCE  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  ANGU. 
119HLAR  POSITION       =»F10.3) 

109  FORMAT(47H  MOMENT  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  ANGU» 
119HLAR  ACCELERATION   *»F10.3) 

110  FORMAT(47H  MOMENT  DERlVrTlVE  FOR  A  UNIT  INCREMENT  OF  ANGU. 
119HLAR  VELOCITY       =.F10.3) 

111  F0RMAT(47H  MOMENT  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  ANGU, 
119HLAR  POSITION       =.F10.3) 

112  F0RMAT(47H  MOMENT  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  VERT. 
119HICAL  ACCELERATION  =.F10.3) 

113  FORMATU7H  MOMENT  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  VERT, 
119HICAL  VFLOCITY      =,F10.3) 

114  FORMAT(47H  MOMENT  DERIVITIVE  FOR  A  UNIT  INCREMENT  OF  VERT. 
119HICAL  POSITION      =.Fl0.3) 

115  FORMAT(22H  STABILITY  INDICATOR  =,E15.5////) 

116  FORMAT  (  If 'H  ESTIMATED  EHP  =.F20.2///) 
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CALL    EXIT 
END 
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SUBROUTINE  ANGLES 

DIMENSION  TK15)  tVALUE(l5) 

COMMON  AA»  Allt  A22»  ASP*  BB»  lill.  B22.  BEAM*  BE7*.  CC  • 
lCllt  C22,  CLB.  CLO,  C6»  Cv ,  DD,  Dili  D22.  DlFBli  DIFCli 
2DRAG»  EE*  Ellt  F22.  EPSIL,  Ff  Gil*  G22.  LIST1*  LIST2*  NRt 
3RHO*  S»  TAU*  TRIM*  Tl,  VALUE*  V,  VCG*  VKTS.  VM.  W»  WETL* 
4WCHINE.  WKEEL*  YI 

CV=V/SQRTF( 3  2.2*BEAM) 

CLB=2.*W/(RHO*<  (V*BEAM)**2  )  ) 

BETAI=BETA*57.2956 

ICLO=l 

5  CLOl=ICLO 
CLO=CLOl*0.01 

IF (CLO-0.006  5*BETAI*CLO**.6-CLB) 2*3.4 

2  ICLO=ICLO+l 
GO  TO  5 

4  CLO=CLO-.001 

IF(CLO-0.006  5*BETAI*CLO**.6-CLB)3»3.6 

6  GO  TO  4 

3  CLO=CLO 
N=l 

NR  =  N 
Tl (N)=l. 

IF(LIST2-1 ) 11  .1.11 
1  PRINT  1000, CLB. CLO. TKN)  .N 
11  CALL  FACTOR 

IF(ABSF(VALUE(N) ) - . 0001 ) 10  .  10  .  7 
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7  N  =  2 

NR  =  N 
Tl (N)=2.0 
IF(LIST2-1 )12»13»12 

13  PRINT  1000.CLB.CLO.THN)  iN 
12  CALL  FACTOR 

IF(ABSF(\ ^LUE(N) J -.0001 ) 10 t 10 t 8 

8  CONTINUE 

DO  9  N  =  3  ,8 

NR  =  N 

SLOPE  =  (VALUE(N-l)  -  VALUE ( N-2 ))/< Tl ( N-l  )  -TKN-2)) 

TKN)  =  TKN-1)  -  VALUE(N-l  )/SLOPE 

IF(LIST2-1 ) 15.14,15 

14  PRINT  1000.CLB.CLO.TKN)  »N 

15  CALL  FACTOR 
IF(ABSF(VALUE(N))-.0001)10»10,9 

9  CONTINUE 

10  TERM=.5*SINF(BETA)*COSF(TAU)/3.1416/SINF(TAU) /COSF(BETA) 
WKEEL=BEAM*< ASP+TERM) 
WCH I NE=BEAM*( ASP-TERM) 
DRAFT=WKEEL*SINF(TAU)     . 
WETL=ASP*BEAM 
TRIM=TAU*57.2956 
S=WETL*BEAM/COSF(BETA) 
PRINT  1002 

PRINT  1003. V, TR IM, ASP ,WKEEL* WCH I NE»S, DRAG* DRAFT,  WETL 
1000  FORMAT(26H  CALLING  FACTOR  WITH  CLB  =,F9.4,3X,7H  ,CLO  =. 
1F9.4,3X,13H  .AND   TRIM  = » F20. 4 , 15X  ,4H  N  =  .13) 
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1002  F0RMATO1H  EQUILIBRIUM  PLANtNG  CONDITIONS//) 

1003  F0RMATU7H  VELOCITY  (FPS)  =.F7.2/14H  TRIM  (DEG.)  =»F6.2/ 
H5H  ASPECT  RATIO  =,F5.2/19H  WETTED  KEEL  (FT)  =»F6.2»20X» 
220H  WETTED  CHINE  (FT)  =,F6.2/22H  WETTED  AREA  (FT**2)  «. 
3F7.2.12X»13H  DRAG  (LBS)  =»F7.2/13H  DRAFT  (FT)  =»F5.2/ 
426H  MEAN  WETTED  LENGTH  (FT)  =,F6.2///) 

RETURN 
END 
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SUBROUTINE  FACTOR 


DIMENSION  TK15).  VALUEU5) 

COMMON  AA.  All,  A22.  ASP »  BB.  Bill  B22»  BEAM.  BETA.  CC » 
lClli  C22.  CLB.  CLO.  CG,  CV ,  DD,  Dili  D22i  DlFBli  DIFCli 
2DRAG.  EEi  Ell.  E22.  EPSlL,  F.  Gil.  G22.  LISTli  LIST2i  NRi 
3RHOi  S.  TAU.  TRIM.  Tit  VAl_U"i  V.  tfCGi  VKTS.  VMt  Wt  WETLi 
4WCHINE.  WKEEL.  YI 

N  a  NR 

IF  (Tl (N)  )2»2.3 

2  IF(LIST2-1 )22.4.22 
4  PRINT  100.N.T1 (N) 

22  Tl  (TO  =  .5 

3  ASP  =  80.*(CLO/Tl (N)**l.l) 
TAU=T1 (N)/57«2956 

IF( ( .0l2*SQRTF(ASP)+.0055*ASP**2.5/CV**2)*Tl ( N ) **1 • 1-CLO ) 
111  .8,12 

11  ASP=ASP+.0l 

IF( ( .012* SQRTF(ASP)+.0055*ASP*«2.5/CV**2)*T1 ( N ) **1 . 1-CLO ) 
111.8.8 

12  ASP  =  ASP  -  .01 

IF(  (  .012*SQRTF(ASP)+.0055*ASP***2.5/CV**2  )*T1  (  N)  **  1  •  1-CLO  ) 
18.8,12 
8  C«CG- (.75-1 •/( 5*2 1* ( CV/ASP )«*2+2 .39) )*ASP*BEAM 

VM  =  V*SQRTF(1.-.012*TAU**1.1/SQRTF(ASP)*(85.-50«*8ETA)/ 
1C0SF(BETA)**2/C06F( TAU) ) 

RE=131770.*ASP*BEAM*V 

REE=.43429448*L0GF(RE) 
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CF=. 008179 

DEL=.001 
50  CF=CF-DEL 

FRE=.242/SQRTF ( CF ) -. 432944  8*LOGF ( CF ) 

IF(REE-FPE) 55.95,50 
55  IF(DEL-. 001 )65. 60.60 
60  CF=CF+DEL 

DEL=.0001 

GO  TO  50 
65  IF(DEL-. 0001)75. 70.70 
70  CF=CF+DEL 

DEL=. 00001 

GO  TO  50 
75  IF(DEL-. 00001 )85.80»80 
80  CFaCF+DEL 

DEL=. 000001 

GO  TO  50 
85  IF(DEt-. 000001 )95. 90.90 
90  CF=CF+DEL 

DEL=. 0000001 

GO  TO  50 
95  CFT=CF+.0004 

A=VCG-.25*BEAM*SINF(BETA)/C0SF(BETA) 

DRAG  =RH0*( VM*BEAM)**2*ASP*CFT*.5/C0SF(BETA) /COSF(TAU)  + 
1W*SINF(T/  J)/C0SF(TAU) 

VALUE(N)=W*(C*(COSF(EPSlL)-SINF(TAU)*SINF(TAU+EPSIL))-F* 
1SINF(TAU)*C0SF(TAU) )+DRAG*(C*(SINF( TaU)*COSF(EPSIL)- 
2SINF(TAU  +  EPSIL)  )  +  ( A*C0SF ( EPS  I L ) -F ) *C0SF ( TAU )  ) 
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IF (LIST2-1 )6»5  .6 

5  PRINT  101  .VALUE(N)  »CFT 

100  F0RMAK47H  NEGATIVE  ANGLE  ENCOUNTERED  ON  ITERATION  NUMBERt 
1I3»40H  THE  VALUE  OF  THIS  ANGLE  (IN  DEGREES)  IS»F10»3) 

101  F0RMATO2H  RETURNING  TO  ANGLES  WITH  VALUE- ♦ E12 . 5 ♦ 
110H  AND  CFT  =,E12.5///) 

6  RETURN 
END 
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SUBROUTINE  COEFF1 

DIMENSION  TK15).  VALUEU5) 

COMMON  AA.  All.  A22.  ASP»  BBt  Bll»  B22»  BEAM»  BETA»  CC» 
1C11»  C22.  CLB.  ClO.  CG»  CV ,  DD ♦  Dili  D22i  DIFB1.  DIFCli 
2DRAG»  EE.  Elli  E22.  EPSIL,  Fi  Glli  G22.  LIST1.  LIST2»  NR. 
3RHO»  S»  TAU.  TRIM.  Tit  VALUE.  V.  VCG,  VKTS.  VM»  W.  WETL. 
4WCHINE.  WKEEL.  YI 
C 

C      EQUATIONS  OF  MOTION  OF  PLANING  HULLS 
C      AA.BB.ETC=ROUTH  CRITERION  FACTORS 

C      A11»B11.ETC=N0NDIMENSI0NAL  FORCE  AND  MOMENT  COEFFICIENTS 
C      Al »B1 .ETC=FORCE  AND  MOMENT  COEFFICIENTS 
'C 

EPS  =  EPSIL 
A  =  ASP 

VERM=.125*RHO*3.1416*(WETl**2)*6EAM 
VERYI=.062  5*.062  5*RHO*3.1416*< WETL**4 ) *BEAM 
CPL = WET L*( 0.75-1.0/ ( 5.21*( ( CV/A ) **Z ) +2 . 39 ) ) 
C 

C      DIFFERENTIATION  OF  LIFT  COEFFICIENT  X  AREA  WITH  RESPECT 
C      TO  TAU  - 


C 


TAU1=TAU-0.001 
TAU2=TAU+0.001 
WC  =  WCHINE 
WK  =  WKEEL 
A  =  l./ASP 
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C 

c 
c 


CLV0L1=1./(2.*WETL*(CV**2 ) )*( ( ( WC**2 ) *S I NF ( 2.*TAU1 ) 
1/BEAM+1./3.*(2.*WC+WK)*SINF(BETA)/C0SF(BETA) ) ) 
3  CLB1=0.5*CLVOL1 

CLSA1=( I4 :570  8*A*TAU1*(C0SF t TAU1 )**2)*( 1  «0-S I NF ( BETA ) ) 
1/(1.0+A)+4.0*(SINF(TAU1)**2)*(COSF(TAU1 )**3)*COSF(BETA)/ 
23.0  +  C|_Bl)*S 

CLVOL2=l./(2.*WETL*(CV**2) )*(  ( < WC**2 ) *S I NF ( 2 • *TAL 2 ) 
1 /BEAM+1 • /3. *( 2.*WC+WK ) *S INF (BETA) /COS F( BETA) ) ) 
5  CLB2=0.5*CLVOL2 

CLSA2=( 1.5  708*A*TAU2*(COSF(TAU2)**Z)*(1«0-SINF(BETA) ) 
l/( 1 .0  +  A) +4.0* (SINF(TAU2)**2)*(COSF(TAU2)**3)*COSF( BETA)/ 
23.0+CLB2)*S 

CL=1.5  708*TAU*A*COSF( TAU)**2*( 1 .-S INF ( BETA) )/(l«+A)+4»* 
1SINF(TAU)**2*C0SF( TAU ) **3*COSF ( BETA ) /3.+ ( CLVOL1+CLVOL2 ) /4, 

TA0=TAU1 

TAB=TAU2 

CO=1.5708*TAU*A*COSF( TAO)**2*( 1 .-SI NF ( BETA ) )/(l.+A>+4«* 
1SINF(TA0)**2*C0SF(TA0)**3*C0SF(BETA) /3.+(CLVOLl)/2. 

CB=1.570  8*TAU*A*COSF( TAB ) **2* ( 1 . -S I NF ( BETA ) )/ ( l.+A)+4.* 
1SINF( TAB )**2*COSF( TAB )**3*COSF( BETA )/3.+(CLVOL2)/2. 

DADTAU=1000.*W/(RHO*V**2 ) * ( 1 . /CB-1 . /CO ) 

DIFB1  =  (.CLSA2-CLSA1  )  /  •  002+CL*DADT  AU 

DIFFERENTIATION  OF  LIFT  COEFF  X  AREA  WITH  RESPECT  TO  Z 

DELZ=0.001 

DELA=-DELZ*BEAM/SINF(TAU)/WETL**2 
DELS  =  BEAM*DEl_Z/(SINF(TAU)*COSF(BETA) ) 
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DELL=DFLZ/SINF(TAU) 

CLV0L3=( ( (WC-DELL)**2)*SINF{TAU) /BCAM+ ( 2 ,*WC+WK-3 .*DELL) 
1*5  INF (BETA)/ ( COSF ( BET A ) *3 . )  ) / ( 2 . * ( WETL-DELL) * ( CV**2 ) ) 
7  CLB3=0.5*CLVOL3 

CLSA3=( 1.5708*(A-DELA)*TAU*(COSF(TAU) **2)*( 1 • -S I NF ( BETA ) ) 
1/ (l.+A-DELA  )+<*.*(  SI  NF(TAU)**2  )*(  COS  F  (  TAU  )  **3  )  *COSF  (  BETA  )  /3  • 
2+CLB3)*(S-DELS) 

CLV0L4=( ( (WC+DELL)**2)*SINF(TAU)/BEAM+(2.*WC+WK+3.*DELL) 
1*SINF(BETA) /( COSF ( BETA )*3. )  )  / ( 2 . * < WETL  +  DELL) * ( CV**2 ) ) 
9  CLB4=0.5<  CLV0L4 

CLSA4=<  1«?708*(  A  +  DELA)*TAU*(COSF(  TAU)**2  )*( 1 .-S INF < BETA ) ) 
l/( 1.+A+DELA)+4.*(SINF(TAU) **2)*(COSF( TAU ) **3 ) *COSF ( BETA ) /3 • 
2+CLB4)*(S+DELS) 

DIFC1=(CLSA4-CLSA3) /0.002 

DIFFERENTIATION  OF  MOMENT  WITH  RESPECT  TO  Z 

CPL1=. 7 5- (WETL-DELL) /( 5 . 2 1 *( CV*BEAM ) **2 /( WETL-DELL ) **2+2 . 39 ) 

CPL2=.7  5-(WETL+DELL)/(5.2l*(CV*BEAM)**2/ (WETL+D ELL) **2+2. 39) 

DCPLDZ=(CPL2-CPLl)/.00  2 

C1=0.5*RHO*( V**2)*DIFC1 

DMDZ  =(  C.G-C PL  )  *C1 /COSF  (  TAU  )-W*DCPLDZ /COSF  (TAU) 

Al=W/32.2+VERM 

B1=0.5*RHO*V*DIFB1 

C1=0.5*RHO*( V**2)*DIFC1 

D1=VERM*VCG* (COSF(TAU)*(CG-.5*WETL) /VCG-SINF( TAU) ) 

El=Bl*(COSF(TAU)*(CG-WETL/2* ) -VCG*S INF ( TAU ) ) 

G1=V*B1 
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A2=YI+VERYI 

B2=E1*( (CG-CPL)*COSF(TAU)-VCG*SINF(TAU) ) 

C2=G1*< (CG-CPL)*COSF(TAU)-VCG*SINF( TaU) ) +W* ( CPL* ( S INF ( TAU ) -VCG* 
1C0SF(TAU)/SINF(TAU)/WK*(C0S-(TAU)/SINF(TAU)-1./SINF(TAU) ) )-CG* 
2SINF( TAU)-VCG*COSF(TAU) ) 

D2-aVERM*(COSF(TAU)*(CG-WETL/2.  ) -VCG*S INF ( TAU )  ) 

E2=B1*( (CG-CPL)*C0SF(TAU)-VCG*SINF(TAU) ) 

G2  =  DMDZ 

E1=E1+V*VERM 

G1=G1+V*B1 

B2=B2+V*D2 

C2*C2+V*E? 

U  IS  AN  ARBITRARY  NOND I  MENS  I ONAL I Z I NG  VELOCITY 
U  =  10. 

QUANTITY  0.5  RHO  CANCELLED  OUT  OF  ALL  FOLLOWING 


A11=A1/(3EAM**3) 

B11=B1/(U*(BEAM**2 ) ) 

C11=C1/(BEAM*(U**2) ) 

DllvDl/(BEAM**4) 

Ell=El/(U*(BEAM**3) ) 

Gll=Gl/( (3EAM*U)**2 ) 

A22=A2/(8EAM**5) 

B22=B2/(U*(8EAM**4) ) 

C2  2=C2/( (3EAM**3)*(U**2) ) 

D22=D2/(8EAM**4) 
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E22=E2/(U*(BEAM**3) ) 
G22=G2/( <BEAM*U)**2 ) 
AA=1. 

BB=( A22*B11+A11*B2  2-D22*E11-E22*D11)/(A11*A2  2-D1 1*022) 
CC*(A22*C11+B22*B11+A11*C22-D22*G11-E22*E11-G22*D11)/(A11 
1*A22-D11*D22) 
DD=(B2  2*C11+B11*C2  2-E2  2*G11-G22*E11  )  /  (  Al  1*A2  2~D1 1*022  ) 
EE=(C2  2*C11-G2  2*G11 ) / ( Al 1*A22-D1 1*D22 ) 
RETURN 
END 
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